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Abstract 

The review is based on the author's papers since 1985 in which a new approach to the 
separation of variables (SoV) has being developed. It is argued that SoV, understood gen- 
erally enough, could be the most universal tool to solve integrable models of the classical 
and quantum mechanics. It is shown that the standard construction of the action-angle 
variables from the poles of the Baker-Akhiezer function can be interpreted as a variant 
of SoV, and moreover, for many particular models it has a direct quantum counterpart. 
The list of the models discussed includes XXX and XYZ magnets, Gaudin model. Nonlin- 
ear Schrodinger equation, 5'L(3)-invariant magnetic chain. New results for the 3-particle 
quantum Calogero-Moser system are reported. 
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1. Introduction 

The separation of variables (SoV), at least, in its most elementary forms such as 
SoV in cartesian, spherical or ellipsoidal coordinates, is an indispensable part of the basic 
mathematical/physical curriculum. Briefiy, SoV can be characterized as a reduction of a 
multidimensional problem to a set of one-dimensional ones. Originated from the works of 
D'Alembert and Fourier (wave theory) and Jacobi (Hamiltonian mechanics), the SoV for 
the long time was the only known method of "exact" solution of problems of mathematical 
physics. However, in the last decades the new techniques, including Inverse Scattering 
Method (ISM) as well as its quantum version (QISM) together with Bethe Ansatz, seemed 
to oust the SoV out of fashion. 

The aim of the present review is to draw attention to the recent progress in under- 
standing SoV and its relations to ISM and QISM. I am going to argue that SoV is far yet 
from being outdated and, even more, has good chances to remain as the most universal 
method of solving completely integrable (classical and quantum) models. There are two 
basic observations which give support to that claim. First: for the classical integrable 
systems subject to ISM the standard construction of the action-angle variables using the 
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poles of the Baker- Akhiezer function is in fact equivalent to a separation of variables. And 
second: in many cases it is possible to find the precise quantum analog of this construction. 

This point of view has being gradually clarified since my first publications of 
1985 which were deeply infiuenced by I. V. Komarov (St. Petersburg State University) 
who used to insist that SoV and ISM are closely related and to underline that the notion 
of SoV should not be restricted to the purely coordinate change of variables as it was 
frequently done. Generally speaking, SoV can be produced by a complicated canonical 
transformation involving both coordinates and momenta. In the quantum case the canon- 
ical transformation should be replaced with a unitary operator. The papers of Gutzwiller 
on 3,4-particle Toda lattice and of Komarov on Goryachev-Chaplygin top Q had pro- 
vided the first examples of how such a unitary operator could be guessed using known 
classical SoV. In |]T|J^ the results of Komarov and Gutzwiller were reproduced using the 
machinery of QISM 0. The algebraic techniques of QISM (R-matrix method) opened the 
way to methodical construction of quantum SoV for the whole classes of integrable models 
generated by different R-matrices (solutions to the Yang-Baxter equation). Though the 
realization of this program is far yet from completion, the list of the models, to which the 
new approach has been applied successfully, grows steadily and presently includes: 

• XXX magnetic chain [^,0] 

• XXX Gaudin model ||], see also [|^ 

• XYZ magnetic chain (classical case) [|T7| 



Goryachev-Chaplygin top 0,0, see also 



Toda lattice g , including relativistic case and boundary conditions |2y] ; see also 
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• Nonlinear Schrodinger equation (infinite volume) ||2z 

• sinh-Gordon model (infinite volume) |j2^ 

• SL{3) magnetic chain |2|,|5l 

• 3-particle Calogero-Moser model, see Section 7 of the present paper 

The early stages of the work were summarized in the reviews [^0], where the term 
"Functional Bet he Ansatz" was used instead SoV, the latter one seeming now to be the 
more accurate. I tried to avoid in the present review the detailed discussions of particular 
models which one can find in but, instead, to summarize the development of the field 
using the unifying concept of the Baker-Akhiezer function. The importance of choosing a 
proper normalization of B-A function is stressed. More attention is paid to the particular 
models which were not discussed in the previous reviews (Nonlinear Schrodinger equation 
in the infinite volume, classical XYZ magnet). New results are reported concerning SoV 
in the quantum 3-particle Calogero-Moser model (elliptic and trigonometric). 
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2. SoV: general notions 

2.1. Basic definitions 

Let us start with tlie classical case. Consider a Hamiltonian mechanical system having 
a finite number D of degrees of freedom and integrable in the Liouville's sense. It means 
that one is given a 2L)-dimensional symplectic manifold (phase space) and D independent 
hamiltonians Hj commuting with respect to the Poisson bracket 

{if„ifj = j,k = l,...,D. (2.1) 

A system of canonical variables (p, x) 

{xj.Xk} = {pj.Pk} = {Pj,Xk} = Sjk (2.2) 

will be called separated if there exist D relations of the form 

^J{xJ,pJ,H^, H2,...,Hd) = i = 1, 2, . . . , (2.3) 

binding together each pair (pj , Xj ) and the hamiltonians . Note that fixing the values 
of hamiltonians = one obtains from ( p.3|) an explicit factorization of the Liouville 
tori into the one-dimensional ovals given by equations 

$j(xj,Pj,/ii,/i2,...,/iD) = j = l,2,...,D. (2.4) 



The action function 5'(/i, x) is defined as the generating function of the canonical 
transformation from (p, x) to the action-angle variables (/, 0) that is 

a.Hii), P,-'-^^, ^^^esiHM. (2.5, 

and satisfies the Hamilton-Jacobi equation 

H,\?^^..-]^k,H) (2.6) 



for each Hj considered as a function of (p, x) . The relations (|2.4| ) allow immediately to 
split the complete solution to partial differential equation ( |2.6|) into the sum of terms 



S{h,xi, . . . ,xd) = Si(h,xi) -\ \- Snih^xn) (2.7) 
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satisfying each the ordinary differential equation 



$,(x„ d^^S,, hi, h2,...,hD) = (2.8) 

(we consider h in Sj{xj) = Sj{h, fixed parameters). 

The last result justifies apparently the term "separation of variables" . It is important 
to warn the reader that many authors restrict the notion of SoV only with the situations 
when the separated variables (p, x) are obtained from some original canonical variables, 
say (P, Q), by purely coordinate transform 

x,=x,m pj=J2^p^- (2-9) 

k 

That condition, though reasonable for a certain class of problems, would be, however, 
too restrictive for our purposes excluding almost all the examples we are going to consider. 

Let us examine now the quantum case. The condition (|2.1| ) is replaced by the com- 
mutativity of quantum operators 

[Hj,Hk]^0 j,k^l,...,D (2.10) 

the relations ( |2.3| ) retaining the same form 

^jixj,Pj,Hi,H2,...,HD)=0 (2.11) 

with the only difference that now we have to fix some ordering of the noncommutative 
operators Xj, pj, H^- Let us assume that the operators in ( p.ll|) are ordered exactly as 
they are enlisted that is x always precedes p and p precedes H (the relative ordering of 
Hn is of no importance since they are commutative ( |2.10| )). 



It is convenient to work in the x-representation that is to realize the quantum states 
as the functions '^{x) from some Hilbert space Ti. The canonical operators Xj and pj 

[xj.Xk] = \pj,Pk] = \Pj,Xk] = -ihSjk (2-12) 

can be realized respectively as the multiplication and differentiation = —ihdj;. operators. 
Let \1/ be a common eigenfunction 

ifj* = hj'i' (2.13) 



of the commuting hamiltonians. Then, under assumptions made about the operator or- 
dering and the reahzation of (p, x) it foUows from ( |2.11j ) that ^' {x) satisfies the equations 



^j(xj, -ihd^^.hi, /i2, . . . , ho)^ = (2.14) 
which suggests immediately the factorization of \I/ 

D 

^{xi,...,xd) = Yl^j{xj) (2.15) 

into functions t/jj [xj ) of one variable only satisfying the equation 

^j{xj, -iWj:.,hi, /i2, • • • , hD)tpj{xj) = 0. (2.16) 

In complete analogy with the classical case, the original multidimensional spectral 
problem (|2.13| ) is reduced to the set of the one-dimensional multiparameter spectral prob- 
lems ( p.l6|) . For determining the admissible values of D spectral parameters hj one has 
thus the system of D equations ( |2.16|) supplemented with appropriate boundary conditions 



determined by the Hilbert space Ti. The function $j in (|2.14| ) can be thought of as a sym- 
bol of a pseudo differential operator. In particular, (|2.16| ) becomes an ordinary differential 



equation if $j is a polynomial in p or a finite- difference equation if $j is a trigonometric 
polynomial in p. 
Semiclassicaly, 

ijf (f ) = e ^ -^(^^ , ipj (xj ) = e^^^ ^""^ ) 
which provides the correspondence of the formulas (|2.7| ) and ( |2.15|) . 

2.2. Magic recipe 

Let us return now to the classical case and discuss the question how to find a SoV for 
a given integrable system. In the XlXth and the beginning of the present century for a 
number of models of classical mechanics, such as Neumann model or various cases of the 
rigid body motion, the SoV was found by guess or some more or less ad hoc methods. 

Here we shall discuss the latest and seemingly the most powerful method based on the 
Baker-Akhiezer function. Suppose, as it is always done in the Inverse Scattering Method 
(ISM), that our commutative hamiltonians Hj can be obtained as the spectral invariants 
of some matrix L{u) of dimensions NxN, called L or Lax operator, whose elements are 
functions on the phase space and depend also on an additional parameter u called spectral 
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parameter. It means that Hj can be expressed in terms of the coefficients tn{u) of the 
characteristic polynomial W{z, u) of the matrix L{u) 

N 

W{z, u) = det{z - L{u)) = 

n=0 (2.17) 

n 

to{u) = 1, tn{u) = tr ^L{u), t N (u) = det L{u). 
The characteristic equation 

W{z,u) = (2.18) 

defines the eigenvalue z{u) of L{u) as a function on the corresponding A'^-sheeted Rieman- 
nian surface of parameter u. The Baker- Akhiezer function Q{u) is defined then as the 
eigenvector of L{u) 

L{u)n{u) = z{u)n{u) (2.19) 

corresponding to the eigenvalue z{u). 

Since an eigenvector is defined up to a scalar factor, to exclude the ambiguity in the 
definition of 0(w) one has to fix a normalization of 0(w) imposing a linear constraint 

N 

Y,0^niu)nM = l (2.20) 

n=l 

where cxniu) may, generally speaking, depend also on the dynamical variables. A nor- 
malization being fixed, 0(w) becomes a meromorphic function on the Riemannian surface 



The poles xj of the Baker- Akhiezer function play an important role in ISM [^. In 
particular, the time evolution of Xj for the hamiltonian fiow with any of the commuting 
hamiltonians can be expressed explicitely in terms of the Riemannian theta-functions 
corresponding to the spectral curve ( |2.1(j| ). Moreover, it was observed that for many 
models the variables Xj Poisson commute and, together with the corresponding eigenvalues 
Zj = z{xj) of L{xj), or some functions pj of Zj, provide a set of separated canonical 
variables for the Hamiltonians Hn- It should be mentioned that, though the seminal 



papers contain all the necessary results, the possibility of their interpretation in terms 
of SoV was not recognized at that time. 

One of the reasons why the poles of Q{u) could provide a SoV is easy to understand. 
Since Zj = z{xj) is an eigenvalue of L{xj) the pair {zj, Xj) should lie on the spectral curve 

W{zj,xj) = 0. (2.21) 
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It remains to observe that, if Zj is a function of pj, the equation ( 2.21 ) fits exactly 
the form ( |2.3| ) since the coefficients tn{xj) of the characteristic polynomial ( |2.17| ) contain 
nothing except Xj and the hamiltonians Hk- 

However, the relations ( |2.3| ) alone are not enough to produce SoV. It is necessary 
that, in addition, the number of the poles Xj be exactly the number of degrees of freedom 



D and that the variables {pj{zj), Xj) be canonical ( p.2|) . Those properties are by no means 
obvious and the last one is rather difficult to verify. To perform the calculation, it is 
necessary to transform the above definitions of {zj,Xj) into a more convenient form. Let 
= res 0(u). From ( |2.19| ) and ( p.20|) there follow, respectively, the eigenvalue equation 

U = Xj 

and the normalization condition for Q^^^ 

N 

L{xj)n^^^ = z^^^^\ a„(xj)Oi^) = 0. (2.22) 



n=l 



Let d{u) denote the one-row matrix [ai{u), . . . ,aN{u)). The existence of a nonzero 
solution O^-') to the problem ( |2.22|) is equivalent to the condition 



rank 



a{Xj) 
L{xj) - Zjl 



N -1 



(2.23) 



which, in turn, can be expressed generically as vanishing of any of two minors of order N, 
for instance 



ai{x) a2{x) <y3{x) 

L2l(x) L22{x) - Z L23{x) 
L3i{x) L32{x) L33{x) - Z 



Lni{x) Ln2{x) 
ai{x) a2{x) 
Lii{x)-z Li2(x) 

L3i{x) L32{x) 



LNsix) 
a3{x) 
Li3ix) 

L33ix) - Z 



aN{x) 

L2n{x) 
L3n{x) 



(x) 



(2.24a) 



aN{x) 
Lin{x) 

L3n{x) 

Lnn{x) - z 



0. 



(2.24b) 



Lni{x) Ln2{x) Ln3{x) 

The pairs {zj,Xj) are obtained then as the roots of the system of equations (2.24) 
which allows, in principle, to count them and to calculate the Poisson brackets between 
them. 

In the sections devoted to the nonlinear Schrodinger equation and XYZ model we shall 
need the formulas for the N = 2 case. The equations (2.24) become then quite simple: 



ai[x) 
L2i(x) 



a2[x) 
L22(x) - 



0, 



ai{x) 
Lii(x) - z 



a2{x) 
Li2{x) 



0. 



(2.25) 
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Eliminating z one obtains the equation for Xj 



B(xj)=0, B = a'iLi2 - aia2{Lii - L22) - alL2i. (2.26) 
The eigenvalue Zj is obtained then from 

\ ^ / U = Xj \ / tt = Xj 

(for the sake of brevity we have omitted the argument u in L{u) and a(w) in the last two 
formulas) . 

Generally speaking, there is no guarantee that one obtains from (2.24) the canonical 
P.b. ( |2.2| ) for some pj(zj). Amazingly, it turns out to be true for a fairly large class of 
integrable models, though the fundamental reasons responsible for such effectiveness of 
the magic recipe: "Take the poles of the properly normalized Baker- Akhiezer function and 
the corresponding eigenvalues of the Lax operator and you obtain a SoV", are still unclear. 
The key words in the above recipe are "the properly normalized" . The choice of the proper 
normalization d{u) of Q{u) can be quite nontrivial (see below the discussion of the XYZ 



magnet) and for some integrable models the problem remains unsolved ^9 . 



2.3. r -matrix formalism 

Given a particular L operator and a normalization of O, one is able, in principle, to 
calculate from (2.24) the Poisson brackets for [zj, Xj) though it could be a formidable task. 
There are, however, techniques which simplify the calculation and allow to verify SoV for 
whole families of L operators instead of handling them individually. 



According to a remarkable theorem proved by Babelon and Viallet the commu- 
tativity of the spectral invariants tn{u) ( p. 17 ) of the matrix L{u) 



{tmiu),tn{v)} = (2.28) 

is equivalent to the existence of a matrix ri2{u, v) of order A^^ x A^ such that the Poisson 
brackets between the components of L are represented in the commutator form 

{L{u),L{v)}=[ri2{u,v),L{u)]-[r2i{v,u),L{v)] (2.29) 

1 2 

where the standard notation is introduced [0: L = L ® 1, L = 1 L, r2i(w, f) = 
Vri2{u, v)V, and V is the permutation operator: Vx ® y = y ® x, Vx, y G C^. 
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Generally speaking, the matrix ri2(w, v) is a function of dynamical variables. So far, 
very little is known of such r-matrices, apart of few particular examples pT|-|3^ . The best 
studied one is the case of purely numeric (c-number) r-matrices satisfying the classical 
Yang-Baxter equation 

[ri2{ui,U2),ri2,{ui,U'i) +r23(M2,W3)] - K3(W1,W3),?'32(W3,W2)] = (2.30) 



which ensures the Jacobi identity for the Poisson bracket (|2.29|) , and especially, the case 
of unitary numeric r-matrices, satisfying, in addition, the relation 

ri2(wi,W2) = -r2i{u2,ui) (2.31) 



and depending on the difference ui — U2- For such r-matrices the relation (|2.29|) takes the 
form 

{L{u)Mv)] = [ri2{u~v)X{u)+L{v)] (2.32) 

and ( 12.301 ), respectively, 

[ri2(w), ri3(w + v)] + [ri2(w), r23(i')] + [ri2,{u + v),r22,{v)] = 0. (2.33) 

To a unitary numeric r-matrix one can associate not only the Poisson algebra ( |2.32| ) 
whose right hand side is linear in L but also the algebra 

{L{u)X{v)] = [ri2{u - v)Mu)L{v)] (2.34) 

with the quadratic r.h.s. Formally, ( |2.34|) can be put into the form ( |2.29| ) with the dynamic 

r-matrix ri2(w, v) = |(ri2(w - v)L{v) + L{v)ri2{u - v)), r2i{u, v) = -^{ri2{u - v)L{u) + 
1 _ 

L{u)ri2{u — f)), but the very special structure of ri2 allows to consider the formula ( |2.34| ) 
rather as a modification of ( |2.32|) . Obviously, (|2.32|) can be obtained from ( p. 34] ) if one 
substitutes L := 1 -\- eL -\- 0{e^), r := er and let £ — > 0. 

Another example of the quadratic P.b. algebra associated to a unitary numeric r- 
matrix is the algebra [p5 | 

12 121 22 1 

{L{u),L{v)}= [ri2{u-v),L{u)L{v)]+L{u)ri2{u + v)L{v)-L{v)ri2{u + v)L{u). (2.35) 

There exists a profound algebraic theory of the unitary numeric r-matrices which 
allows to classify r-matrices in families labelled by Lie algebras. A particularly important 
example is given for any semisimple Lie algebra Q by the formula 



ri2(t^) = - V/a®/a (2.36) 

II ' ^ 
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where p is a numeric constant and /q, G ^ is an orthonormal basis with respect to the 
KiUing form. The result does not depend on the choice of the basis. Taking then various 
finite-dimensional representations of Q for the generators one can obtain from ( |2.3(j| ) 
the family of r-matrices related to Q. 

In particular, for Q = gl{N) and the fundamental vector representation, one has 

ri2{u) = ^V. (2.37) 

The last example deserves special attention since, so far, it is the only series of r- 
matrices for which a general SoV construction is obtained. 



3. G'L(iV)-type models 
3.1. Classical case. 



It turns out that in case of the (7L(A^)-invariant r-matrix ( |2.37| ) the normalization 
of the Baker- Akhiezer function 0(u) corresponding to any constant numeric vector a. in 
( p.20|) produces SoV. The simplest choice of a is 



aN-i{u) = 0, ctAT = 1 



VLn{u) = 1, 



0. 



(3.1) 



The corresponding separated coordinates Xj are defined as the poles of Q(u), and the 
canonically conjugated momenta ( |2.12| ) are pj = —p~^Zj for the linear P.b. ( p.32|) and 
Pj = —p~^\nzj for the quadratic P.b. (|2.34| ). For the linear P.b. case the above results 
were obtained in The quadratic P.b. case was studied in for = 2, 3 and 

generalised to arbitrary N in |^ . 

In case of the normalization ( |3.1| ) the equations (2.24) for [z, x) simplify a bit 



(3.2a) 



L2i{x) 


L22{x) - Z 


L23{x) 


L2,N- 


-i(^) 


Lslix) 


^32 (a:) 


Ls3{x)-Z . 


L3,N- 




Ln-i,i{x) 


LN-l,2ix) 


LN-l,3ix) . 


■ -^Af,Af-l 


(x) — z 


Lni{x) 


Ln2{x) 


Ln3{x) 


Ln,n 


-i{x) 


Lii(x) - z 


Li2{x) 


Lisix) 






Lsiix) 


L^2{x) 


Ls3{x)-Z . 


L3,N- 


_i(a;) 




Ln-1,2{x) 


LN-l,3ix) . 




(x) — z 


Lni{x) 


Ln2{x) 


Ln3{x) 


Ln,n 


-i{x) 



(3.2b) 
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The equations (3.2) themselves can be used aheady for calculation of P.b. between 
X4 and z J ^ sec [p7| , |38[] . However, they are not convenient for quantization because of the 
operator ordering problem. So, we take one more step and eliminate zj from (3.2). The 
result is one equation 



B(x, 







(3.3) 



for Xj where B{u) is a certain polynomial of degree A^(A^ — l)/2 in components of L{u). 
The corresponding eigenvalue Zj is obtained then as the value 



A{xj 



of certain function A{u) expressed rationally in components of L{u) 
For instance, for N = 2 

B{u) = L2i{u), A{u)=Lii{u). 

Note that for u — xj the matrix L{u) becomes triangular 

which explains why its eigenvalue Zj is given by A[xj). 
For = 3 one obtains 



Lii Li2 
L31 L^2 



B{u) = L^i{u 
There are two possible ways to choose A{u) 

(u) 



u) + Lz2{u) 



A^{u) 



Lii L12 
L31 L32 



-^32 (w) 



A2{U) 



L2I L22 
L31 L32 



L2I L22 
L31 L32 



[u 



which are equivalent modulo B{xj) = 



Ai{xj) = A2{xj) 



Zj = Ai{xj) = A2{xj). 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.^ 



(3.9) 



The expressions for A{u) and B(u) for the general N are given in p3]. We have to 
warn the reader that due to a different choice of normalization of fl{u) the formulas ( |3.5| ), 
]8D differ from those in |l7|, p5| , |38[ . 
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The equations (|3.3| ) and (|3.4|), like (3.2), can also be used for calculating the P.b. 
between Xj and Zj. Their advantage for the sake of quantization is that B{u) Poisson 
commute 

{B{u),B{v)} = (3.10) 

which entrains immediately the commutativity of Xj (see the next subsection). 

The correct P.b. between Xj and zj are not enough to establish SoV. The last condition 
(usually, easy to verify) is the correct number of variables Xj which should be equal to the 
number D of degrees of freedom. In some degenerate cases, the number of Xj could be less 
then D and some additional variables should be added (see example of Calogero-Moser 
model in section 7). 

To conclude this subsection, let us stress that so far no generalization is known of the 
above results to the r-matrices corresponding to the Lie algebras other then the An series. 



The difficulty is that the simplest normalization (|3.1|) does not work more: the function 
B{u) has too many zeroes (more than the number D of degrees of freedom) and they do 
not commute [^. Hopefully, some other normalization ( |2.2UD will work which remains a 



challenging problem. 
3.2. Quantization 

In the quantum case the components Lmn of the NxN matrix L{u) become quantum 
operators, and the Poisson brackets should be replaced by some commutation relations 
sutisfying the correspondence principle [,] = — }. A nice feature of the Poisson 
algebras ( |2.32| ) and (p.34| ) is that (in contrast with the general case ( p.29| ) of dynamical 
r-matrices) it is well known how to quantize them. 

In the linear case ( p.32| ) the quantization is straightforward 

[L{u),L{v)] = [xi2iu-v),L{u) + L{v)], x{u) ^ -ihr{u), (3.11) 
in the quadratic case ( |2.34| ) it is more tricky. The algebra ( |2.34| ) is replaced by 

12 2 1 

Ri2{u - v)L{u)L{v) = L{v)L{u)Ri2{u - v) (3.12) 
where R{u) satisfies the quantum Yang-Baxter equation 

Ri2iu)Ris{u + v)R23{v) = R23{v)Ris{u + v)Ri2{u) (3.13) 
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and is related to r{u) through the semiclassical expansion 



R{u) = l + inr{u) + 0{n^). (3.14) 



In the GL{N) case ( |2.37|) the quantum r and R matrices are 



V V 
t(u) = —r) — , R(u) = l+rj — , rj = iTip. (3.15) 

u u 

The relations ( |3.12|) define the associative algebra Y[gl[Ny\ called yangian of gl{N). 
The quantum integrals of motions tr^iu) for the yangian are obtained by appropriate 
deformation of the classical formulas ( |2.17D 

tn{u) = trL(w) A L{u + 77) A ... A L{u + (n - 1)?]) (3.16) 

[t^{u)Mu)]=^ (3.17) 
The quantity tN{u) {the quantum determinant of L{u)) produces central elements 

[tN{u),L{v)]=Q (3.18) 

of the yangian. Naturally, on the irreducible representations tN{u) is a number- valued 
function. 

The quantum analog of the construction of the SoV based on the equations and 
( p.4|) is found presently only for N = 2 and = 3. The formulas given below are taken 
from [)7|,[23[1 up to small variations due to the different choice of normalization of Q{u). 
In the GL[2) case the operator-valued functions A{u) and B{u) are defined by the same 
formulas (|3.5|) as in the classical case. By virtue of (|3.12|) the operator family B{u) turns 
out to be commutative 

[B{u),B{v)] = (3.19) 

which allows to define the operators x„ as the commuting "operator roots" of the equation 
( p73|) (for the mathematical details see 0). To give a sense to the formula ( p^) in the 
quantum case, it is necessary to fix the operator ordering. Assume that x's in (|3.4|) are 
ordered to the left that is 

A{x^) = J2xjAk for A(w) = ^wMfc. (3.20) 

k k 
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Then, using the commutation relations ( 3.12 



) it is possible to verify the relations 



J' 



Xk] = [zj,Zk] = 0, ZjXk = {xk + r]^jk)zj (3.21) 



which suggest the realization of the operators Zj as the shift operators in an appropriate 
Hilbert space 7i 3 \l/(x) of functions on the common spectrum of the operators xj: 

^,^'(x) =C(x,)^'(...,x,+ 77,...). (3.22) 

The choice of the function ({x) in (|3.22[) is dictated by the properties of the Hilbert 
space TC depending on the concrete model (see examples below). Note that there is certain 
liberty in choosing Ci^j) due to the canonical transformations 

^{x)^l[u{x,)^{x) ^ axj)^'^^^^j^ax,). (3.23) 

The SoV follows then from the relations 

- Zjhixj) +t2ixj) = (3.24) 



which generalize the classical characteristic equation (|2.17 ) and fit the form ( |2.11j ) required 



for the quantum SoV (note the operator ordering!). Denoting by rn{u) the eigenvalues of 
the commuting operators tn{u) one obtains for the corresponding separated equation ( |2.16| ) 
the finite- difference equation of order 2: 

Cixj + r])C{xj)'i/;j{xj + 2ri) - C{xj)Ti{xj + r])'i/;j{xj + ??) + T2{xj)'iljj{xj) = 

or in more symmetric form 

Ti{xj)tpj{xj) = A_{xj)tpj{xj -rj) + A+{xj)tpj{xj + r}) = (3.25) 

where 

A,(u) = CM, A_(«) = f^. 

C{u - r]) 

In the GL{3) case the quantum B{u) and A{u) are obtained as deformations of the 
classical formulas ( |3.7| ) and (|3.8|): 

B{u) = Lsi{u - r])[Ls2iu)Lii{u - 77) - Lsi{u)Li2iu - 77)] 

(3.26) 

+ Ls2{u - 'r])[Ls2{u)L2i{u - r]) - Lsi{u)L22{u - r])], 
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(3.27) 



Ai{u) = L32^(lt)[L32(w)Lii(u -1])- L3i{u)Li2{u - 7]))], 

A2{u) = -L^l{u)[L32{u)L2i{u -T])- Lsi{u)L22{u - 77))], 

Zj = Ai{xj) = A2ixj). (3.28) 

The rest is similar to GL{2) case. The quantum characteristic equation is 

- zh{xj) + Zjt2{xj) - h{xj) = 0, (3.29) 



and the separated equation is now a third-order finite- difference equation . 

There is httle doubt that the above constructions can be generahzed to the arbitrary 
values of A^, though the complicated structure of B{u) prevents rapid progress. 



The analogous results for the linear commutation relations (|3.11|) can be obtained from 



the formulas for the quadratic relations ( ^.121 ) in the limit L := 1 + eL + 0{e ), rj := £77 



e ^ 0. However, though it is clear enough that the expansion in e of the formulas ( |3.16| ) 
for the quantum integrals of motion for the yangian should produce, in principle, some 
commuting hamiltonians for the algebra ( |3.11| ), obtaining explicit formulas for them when 
is arbitrary remains still unsolved problem, to say nothing about expressions for A{u) 
and B{u). In case of the r matrix of the form ( |2.36| ) Feigin and Frenkel has proved for 



any Lie algebra Q that the quantum commuting operators do exist which are deformations 
of the spectral invariants of the classical L matrix though their method of proof does not 
produce any effective formulas. Nevertheless, the simplest integrals of motion are easy to 
produce. Note that ti(zi) = trL(u) is a trivial central element of the algebra ( |3.11| ), so it 
can be safely put to be 0, which corresponds to considering sl{N) instead of gl{N). The 
first nontrivial invariant is quadratic in L{u) and coincides with the classical expression 

i2{u) = U^L'^u), [i2{u),i2{v)] = (3.30) 

(for the general simple Lie algebra Q one should use the corresponding Killing form). 

The above quadratic invariant is enough to serve the sl{2) case. The definitions ( p.5| ) 
of A(u) and B{u) and, respectively, (|3.3| ) of and (|3.4| ) of Xj and Zj remain the same as 
in the quadratic case. The commutation relations ( ^.211 ) and the quantum characteristic 
equation ( p.24|) are replaced, respectively, by 

[xj, Xk] = [zj, Zk] = 0, [zj, Xk] = Zj5jk (3.31) 

and 

z]-i2{xj)=0 (3.32) 
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(here we put r] = 1 for simplicity). Using then the reahzation 



Zj = d^^ + C{xj) (3.33) 
for Zj one obtains for the separated spectral problem the second order differential equation 

i;" + 2Ci;' + {e + C)i^ = T2i;. (3.34) 

It is easy to anticipate that for GL{N) the separated equation should become an A^th 
order differential equation, though the calculation still waits to be done. 

For discussion of SoV for the quadratic algebra ( p.35| ) in the sl{2) case (open Toda 
chain with boundary conditions) see . 

The above scheme of quantum SoV is general enough to serve the variety of quantum 
integrable models obtained by taking concrete representations of the yangian ( |3.12| ) or 



the algebra ( p.ll ). Below we shall illustrate on few examples of GL(2)-type models the 



diversity of possibilities arising. The main problem of adjusting the general scheme to a 
concrete model is to find the spectrum of the commuting operators B{u) and, consequently, 
of Xj and to describe the functional space in which the one- dimensional spectral problem 
(^) or ( |01) has to be solved. 



4. XXX magnetic chain 

The general finite-dimensional irreducible representation of the yangian 3^[s/(2)] is 
realized in the tensor product of D finite-dimensional irreducible representations of the Lie 
algebra sl{2) (classically, D is the number of degrees of freedom) 

(^n)' + 1(3+ S- + S-S+) = Ul^ + 1) 1, 1 2, . . .} ^^'^^ 

and can be written in the factorized form {monodromy matrix) 

Lxxxiu) = KCd{u -5d)... C2{u - 52)C^{u - 5^) (4.2) 

where 

CM=ul + TlSr,, Sn=(^fl f^s), (4-3) 
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and i^T is a constant numeric matrix. It is customary to use the notation T{u) instead of 
L{u) like e.g. in [|^, but here we had to sacrifice it to preserve the coherence of notation. 
The representation given by ( [4.2|) has dimension ^n=l(2^'^ + -'-) parametrized by: 

2x2 matrix K, number D of degrees of freedom, spins and shifts dn- The Casimir 
operator (quantum determinant (|3.16|) ) t2{u) takes the value 

D 

t2 [u) = del K W{u - 5r, - Lv) {u-dn + InV + v) ■ (4-4) 

The corresponding quantum integrable model is called inhomogeneous XXX magnet. 

The parameters of the representation being in the generic position, the representation 
( ^72|) of the yangian turns out to be irreducible, and the spectrum of the operators Xj 
defined from (|3.3[) for B{u) given by (p.5|) turns out to be the finite set M 



spec[xj] = Aj = {Xj ,Xj + ry, . . . , A+ - r], A+}, = 5j ± Ijij. (4.5) 

For the separated finite- difference equations ( |3.25| ) to be well defined on the finite set 
Aj, the coefficients A-|-(a;) must satisfy the boundary conditions 

The most convenient choice of A±{u) is 

D 

A±{u) = K±Y\_{u — Sn T ^nV)^ K+K- — detK. (4.6) 

n=l 

The s/(2) Gaudin model is the degenerate case of XXX magnet obtained in the limit 
1] ^ 0. The corresponding L operator is produced renormalizing L operator (O), putting 
K :— 1 + ry/C, tr/C = 0, and expanding in rj: 

Lxxx{u) 



1 + ?7^Gaudm(w) + 0(77 ). 



The L operator 



D 

LG.udiniu) = /C + V (4.7) 

satisfies the linear commutation relations (|3.11| ) with the r matrix (|3.15| ) for 77 = 1. The 
spectral invariant t2{u) ( p.30|) produces the commuting hamiltonians which are quadratic 
in spin operators 

D D 



n=l n=l 
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D 



m=l 



S„. - dr. 



(4.9) 



The normalization ( [4.6| ) of A±{u) in ( |3.25|) corresponds to the normahzation 

In. o 



C{u) 



D 

E 

n=l 



det/C 



u 



(4.10) 



of C(w) in (13331) , (13341) . 

In the limit r] — > the 2Zj + 1 points of the spectrum Aj of the operator Xj merge 
into one point Xj = Sj of multiphcity 21 j + 1. The space of functions on Aj is understood, 
respectively, as the the ring of polynomials in Xj factorized over the ideal {xj-dj)^^^+^ = 0. 
The spectrum of the hamiltonians ( [1.9| ) is given then by the values of in ( ^.8| ) for which 
the differential equation ( p.34| ) with C(m) given by ( [1.1 0| ) in each of the points x = 
has a regular solution i^j{x) = 1 + J2T=ii'^ ~ ^j)'^''!^^^ see @. 

If one realizes the spin operators 5^ ( |4.1j ) as the differential operators 



— yndy„ Ir 



s, 



+ 



yndy„ 



(4.11) 



then the equation B{xj) = L2i{xj) = defines a "purely coordinate", in the sense (^ 
change of variables {Vn} {^j}- The separated coordinates X j can be described as 
generalized ellipsoidal coordinates, see 0. In fact, all the models allowing SoV in general- 
ized ellipsoidal coordinates, such as Neumann model |jTl|] and its generalizations ||T5|JT4 



or 



Euler-Manakov top [jT^, can be considered as degenerate cases of Gaudin model [|T0|- For 



the generalization of Gaudin model to osp{l\2) Lie superalgebra see [^. An application 
of Gaudin model with L{u) having a second order pole to the atomic physics (Coulomb 



three-body problem) is considered in [|T2 



The particular simplicity of the Gaudin model makes it attractive for rigorous math- 
ematical analysis. Much attention was devoted last years to the study of the connection 
between Gaudin model and Knizhnik-Zamolodchikov (KZ) equations for the correlators in 
conformal field theory. The fact that the eigenvalue problem for the Gaudin hamiltonians 
( [1.9|) coincides with the KZ equations on the critical level can be exploited both to pro- 
duce integral representations for the solutions to KZ equations and a new derivation 
of the formula for the norm of the Gaudin eigenfunctions [^1|]. In the recent paper 
the representation theory for the affine Lie algebras is applied to derive, in particular, 
Bethe equations for the Gaudin model corresponding to arbitrary simple Lie algebra and 
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to reveal thus the algebraic meaning of Bethe ansatz. Hopefully, the methods developed 
will be useful also in understanding the algebraic roots of SoV. 



m 

The XXX model, as well as the s/(2) Gaudin model, presents a convenient possibility 
to compare the results of SoV method with those obtained by the Algebraic Bethe Ansatz 
(ABA) [0. Since the subject is discussed in detail in we present here only the 

summary of the analysis. The SoV and ABA methods lead to the same equation ( p.25|) , 
resp. ( |3.34| ), though its interpretations differ. In SoV method the equation is solved on the 
finite set A = [Jj-^j whereas in ABA is supposed to be a polynomial whose zeroes 
Vm parametrize the Bethe vector 

M 

"^v = n L2i{vm) |0) , |0) = 0. (4.12) 

m=l 

In the SoV approach there is one-to-one correspondence between the solutions of the 
problem (|3.25|) on A and the eigenvectors of the commuting quantum hamiltonians which is 
not the case in the ABA approach where the so-called "completeness problem" arises. The 
SoV method provides the basis for the rigorous analysis of the completeness problem and 
allows to formulate the criterion of completeness. In case of the Gaudin model the criterion 
sounds as follows [§. Let Q{u) be a polynomial solution to the differential equation ( |3.34|) . 
If Q{dj) 7^ Vj = I,...,!) then the corresponding Bethe vector ( [4. 12] ) is nonzero. If, 
however, Q{Sj) = for some j then the corresponding Bethe vector is zero, and the set 
of Bethe eigenvectors is incomplete. Moreover, the nonzero eigenvector corresponding to 
such Q{u) exists if and only if the linearly independent solution of the differential equation 
is regular in the same point x = 5j. 

The power of SoV is revealed most obviously in the cases when the representation of 
the yangian does not possess the highest vector |0) such that L12 |0) = and hence ABA 
cannot be applied. These cases correspond to the infinite dimensional representations of 
sl{2) for the operators (^73|). The corresponding separated wave functions ilj{xj) are 



not more polynomials, and the separated equations (|3.25|) , (|3.34|) should be accompanied 
with some square integrability conditions. Depending on the real form of s/(2) in question 
there is plenty of analytical possibilities. 

For the Goryachev-Chaplygin top the spectrum of xj is real and discrete, and the shift 
rj in ( |3.25| ) is real, see and also fl^ for generalizations. 

In case of the Toda lattice the spectrum of Xj is real and the shifts 77 in (|3.25|) are 
imaginary [|| . See also for the relativistic version and [|2^ for the lattice with boundary 
conditions. 
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A version of noncompact XXX magnet applied recently to describe the QCD in 
the asymptotic high energy regime also does not have the highest weight vector, so the 
SoV is the natural approach to try. 

For the analogous effects in the Neumann model, see [ pT| . 



5. Infinite volume limit 

So far, we discussed only the integrable models with a finite number D of degrees 
of freedom. The passage to D = oo can be made in two ways: either by taking the 
continuum limit or the infinite volume limit. In the continuum limit the representation 
( ^72|) of the monodromy matrix L{u) as a product of local C operators Cniu) is replaced 
by the representation in the form of the ordered exponential [|5| j4^] 

Lf_{u) =:ex$ j^J ^{u,i)di: . (5.1) 

For example, for the nonlinear Schrodinger equation, described in terms of the canon- 
ical fields ^{i), 

[*(0, ^(O] = = 0, [^(0, = - a (5.2) 

acting in the Fock space \E'(^) |0) = 0, the infinitesimal L operator £(u, ^) is given by 

^v^vl/(0 iu/2 J ^^■'^> 

The corresponding quantum monodromy matrix L^^(it) given by ( [5.1| ), where : : 
stands for the normal ordering, satisfies the commutation relations (|3.12|) with the R 
matrix of XXX type ( 3.15 ) and r] = —ic. The quantum determinant t2{u) of L^"*" (m) is 
equal to e~'^^/^, and the trace ti{u) of L^^(u) generates the commuting hamiltonians, in 
particular 

H = J (^'l^'g +c^'***^'*)rf^. (5.4) 



We assume that the fields \I/(^), ^*(0 periodic in ^ with the period ¥ = — 
and the coupling constant c is positive c > (repulsive case). 

It is convenient to choose the following normalization d{u) of the Baker- Akhiezer 



function (ICTD 



ai{u) = 1, a2{u) = i. 
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The reason for such a choice is that the corresponding operator family B{u) (|2.26 ) 



B{u) = -iLii{u) + Li2{u) + L2i{u) + iL22{u) (5.5) 

has the symmetry B{u)* = B{u) and its zeroes Xj are self-adjoint operators. 
It is easy to find that the separated equation ( ^.25| ) takes the form 

Ti{x)ij{x) = e-^^^/2^(x + ic) + e^^^/2^(x - ic). (5.6) 

Since B{u), like L^^'^{u), is not a polynomial but a holomorphic function of u G C, 
it has an infinite discrete set of zeroes xj having asymptotics xj = 2nj /V + 0{j~^) as 
j ^ oo. The necessity to handle the functions of infinite number of variables complicates 
greatly the justification of the standard SoV construction even in the classical case. 

The situation, however, simplifies drastically in the infinite volume limit (we consider 
the zero density case described again by the Fock representation for \E'(^), \E'*(^)). The 
definition of the monodromy matrix L{u) for V = oo needs some caution. If one tries to 



pass to the limit ^± ^ ±00 directly in the expression ( |5.1| ) one finds immediately that 



in order to get a finite expression, it is necessary to cancel the asymptotic behaviour of 



Lf^ (u) introducing the exponential factors 



Lt^iu) = Jhn^e^"<^3^+/2^e+(^)g-...3^-/2^ ^3 ^ ^1 _0^^ ^^^^^ 



which, however, destroy completely the nice commutation relations ( 3.12 ). The solu- 
tion of the problem is to factorize first the finite- volume monodromy matrix L^^ (u) = 
(w)L|'^ (w) and to permute the factors: -^|°(w) = L^° {u)L^^'^{u). This operation does 
not change the quantity ti{u) = trL(w) generating the integrals of motion. The matrix 
L|°(m), however, behaves in the limit ^± ±00 much better: 



where 



1 0\ p = 

r ~ 1 1 



and the upper (lower) sign corresponds, respectively, to Imu > (< 0). 

The scalar factor e^*"^/^ can be cancelled out of L|°(u) since it does not affect the 
relation ( |3.12| ). Hence, we can take for the monodromy matrix in the infinite volume the 
matrix 

Hu) = Li°^P±L+^ . (5.8) 
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The matrix L{u) is analytical in the complex plane of u except the cut along the real 
axis. Its quantum determinant t2{u) is zero. 

It is interesting to compare the above results with those of Kyoto group on the XXZ 



magnetic chain [^. In both cases the infinite volume case is considered, and a sort of 
monodromy matrix L{u) satisfying ( |3.12| ) is constructed. In the XXZ case, however, all 
components of L{u) are integrals of motion whereas in the nonlinear Schrodinger case only 
tr L{u) is one. Moreover, in the XXZ case L{u) is analytic in C whereas in the NLS case 
L{u) has a cut along the real axis. The above distinctions are probably due to the different 
nature of the vacuum state: antiferromagnetic for XXZ and ferromagnetic for NLS which 
in the last case requires introducing the asymptotic exponents e='=*"'^3^±/^. 

Let us consider now how the whole SoV construction is modified in the infinite volume 
limit. As ^± ^ ±00 the zeroes Xj of B{u) accumulate to a continuous distribution 

2n 27r , , 

with the density q{i^)- The Hilbert space H of quantum states is realized then as the space 
of the linear functionals of q(z^), z/ G M which are square integrable 



|VF||2 = j \W[q{iy)]\'^ 5m. 



with respect to the measure 5m which, fortunately, happens to be Gaussian and is char- 
acterized uniquely by the correlator (covariance kernel) 



<9M9M)„ = ^ln(l+(^) 



The representation of B{u) as an infinite product B{u) ~ Hjl^ ~ ^j) replaced by 
the Cauchy integral 

B{u) = ^- exp j - / —q{^) } (5-9) 



where the upper (lower) sign corresponds, respectively, to Im-u > (< 0). 

Though the quantity ti{u) is represented in ?i as a complicated variational operator, 
its eigenfunctions nevertheless can be found exactly and have rather simple structure which 
is natural to consider as a continual analog of SoV. The separated equation ( |5.6| ) is replaced 
respectively by a boundary problem for analytical functions having a cut along the real 
axis. Unfortunately, the complications due to the cut make a brief explanation impossible 



and for the details we refer the reader to the original papers [22 . 



22 



For the nonlinear Schrodinger equation the above SoV procedure can be justified 
rigorously by comparison with the results known from the algebraic Bethe Ansatz ||22|| . 
The analogous construction for the relativistic sinh-Gordon model, though not so rigorous, 
can also be performed and leads to quite reasonable results concerning the spectrum 



of the model. It would be interesting to generalize the results of to the relativistic 
To da field theories. 



6. Classical XYZ magnet 

The XYZ magnet provides an example of quite nontrivial normalization of the Baker- 
Akhiezer function 0(m) necessary to produce SoV. The construction presented below is 



taken from [17|. The model is described in terms of the 2x2 matrix L{u) satisfying the 



relation ( p.32| ) with the r-matrix 

3 

^(^) = "^WaiujaaCTa (6.1) 



1 2 

Wa{u) 

a=l 

where are standard Pauli matrices 



^' = {i J)' ^2 = (° o')' ^^=(o -i)' (^-^^ 

and Wo,{u) are certain elliptic functions whose exact expression is not important for the 
moment. It suffices to remark that r{u) is a meromorphic function on C having simple 
poles on the periodic lattice T = {u E C\u = m + rn; m, n G Z; Imr > 0} and 
possessing the periodicity properties 

112 2 

r(u + 1) = air(u)ai = air(u)ai, 

1 2 2 (6-3) 
r(u + t) = a3r{u)as = asr{u)a3. 

The L operator L{u), in turn, is a holomorphic function characterized by the quasiperi- 
odicity properties 

L{u + 1) = {-l)^aiL{u)ai, L{u + r) = (-l)^e-^"^("+"V3L(^x)cT3 (6.4) 



where D isa, positive integer. The conditions ( |6.4|) determine L{u) up to 4D free parameters 
||46|| which can be considered as the dynamical variables with the Poisson structure defined 
by ( p.ll|) . Since, however, the determinant detL(w) generating the center of the Poisson 
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algebra contains 2D parameters (Casimir functions), their values can be fixed which leaves 
2iI>-dimensional phase space, the constant D being thus the number of degrees of freedom. 

In contrast with the XXX magnet, the normalization d{u) = const of does not 
produce SoV for the XYZ magnet for any a.. The reason is that the corresponding function 
B{u) has no definite quasiperiodicity for the period lattice V and can be characterized only 
by the periodicity properties in 2V 



B{u + 2)= B{u), B{u + 2r) = g-^'^^^^u+Sr)^^ 



u 



from which it follows that B{u) has 4D zeroes in the fundamental region C/2r whereas 
one needs only D separated coordinates Xj. 

The correct normalization []T^ is given by the holomorphic functions d{u) having the 



periodicity properties 

ai{u + 1) = ai{u + r) = -e-^^(^+^-^)/2ai(w), 

(6.5) 

a^iu + 1) = -a2(w), a2iu + t) = e-*^("+^-^)/2^i(w) 

where y is a parameter. The explicit expressions for an{u) can be given in terms of 
theta- functions for the lattice of periods F |T^. The function B{u) corresponding to the 
normalization vector ( |6.5|) is given by the formula ( |2.26| ) and has good periodicity properties 
on the lattice F 

B{u + 1) = {-l)^+^B{u), B{u + t) = (-l)^+ie-'"(^+i)("+")e-^"2/5(^). (g.e) 

Consequently, B{u) is a theta-function of order D + 1 and has D + 1 zeros in the 
fundamental region C/F. It remains to require one superficial zero of B{u) to be a constant 
(c-number) 

B{0 = (6.7) 

which can be considered as the equation determining the parameter y as a function on 
the phase space. The remaining D zeroes of B{u) are candidates for the separated coordi- 
nates Xj. A direct, though cumbersome, calculation of the Poisson brackets [jl^ validates 
successfully the conjecture. 

An open question is the quantization of the above construction. The difficulty is the 



ordering problem in the expression ( |2.26| ) since the quantities an{u) contain dynamical 
variables through their dependence on y. 

The example of XYZ model raises the question if the correct normalization producing 
SoV could be found for other integrable models where the simplest choice d{u) = const is 



known to fail 29 
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7. 3-particle Calogero-Moser model 

In this section we present new results for the classical and quantum 3-particle ellip- 
tic Calogero-Moser model. The model is interesting as an example of SoV in case of a 
dynamical r-matrix. 

In the classical case the model has 3 degrees of freedom and is described in terms of 
the canonical variables (7r^,g^) 

{qm,qn} = {T^m,'n:n} = 0, {tt^ , ?n } = (^mn , m, n = 1, 2, 3. (7.1) 

The commuting hamiltonians are 



H2 = TTi%2 + TTlTTs + 7r27r3 - q'^ {p{qi2) + piQls) + p(?23)) (7.2) 
H3 = TriTTsTTa - fi'^(7rip(Q23) + 7r2p(?13) + 7I"3p(qi2)) 

where Qmn = Qm — Qm p is Weierstrass elliptic function, and g is the coupling constant. 
The corresponding L operator is 

TTi -igQi2{u) -igQisiu) 
^cm(w) = I -igQ2i{u) 712 -wQ22,{u) I , (7.3) 



-igQsliu) -igQ32{u) tts 



where 



a[u)a[qmn) 

and cr{u) is Weierstrass sigma function. 

The hamiltonians (|7.2[) can be obtained from the spectral invariants of the L operator 



det(z - L{u)) = - z^ti{u) + zt2{u) - t3{u) (7.5) 
ti{u) = Hi 

t2iu) = H2+3g^piu) (7.6) 
tsiu) = H3+g^p{u)Hi - ig^p'iu) 
The L operator ( |7.3| ) satisfies the identity ( |2.29[ ) with rather complicated r-matrix 



depending on q [^. In the absence of general theory of dynamical r-matrices the only 
available strategy is to try one-by-one possible ansatze for the normalization d{u) of the 
Baker- Akhiezer function Q{u). Fortunately, the very first attempt succeeds: the simplest 
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normalization (SJ.) which was apphed to the GL(A^) -magnet does also produce SoV for 
the Calogero-Moser model. 

For our purposes it is convenient to write down the set of equations for the pair (x, z) 



as (|3.9|) where the functions Ai^2iu) are given by the formulas ( p?8|) and, in our case, are 
Ai{u) = TTi + ig[C{u) - C{u - Q23) + C(9i2) - Ciqis)] 

(7.7) 

A2{u) =712 + tgiCiu) - C{U - qi3) - C{Q12) " C(?23)] 

where ({u) is Weierstrass zeta function. 

Since the r-matrix for the CM model is different from (|2.37| ) we cannot rely on the 
results obtained for the GL{N) magnet and have to calculate the Poisson brackets between 
z and X directly. It turns out that the equations ( |3.9| ) have only two solutions: {zi,xi) 
and {z2, X2). For the third pair of variables one has to take (P, Q) 

P = ni+n2 + ns, Q = q^. (7.8) 

The calculation shows that the variables (P, zi, Z2; Q, xi, X2) are canonical and satisfy 
the relations 



P-Hi=0 

(7.9) 

4 - z]Hi + Zj{H2 + 3g^p{xj)) - {H3 + g^p{xj)Hi - ig^p'{xj)) = 



which fit the form ( |2.3| ) and provide thus a SoV. 

In the quantum case the momenta ttj are realized as the differentiations TVj = —idq.. 
The hamiltonians (|7.2| ), respectively, are replaced by the differential operators 

H2 = -dlq^ - dlq^ - dlq^ - gig - l)[p((/i2) + P(?i3) + P(?23)], (7.10) 

^3 = ^5^19293 + ^^(f - l)[p(923)9g, + p(?13)9g2 + p{(l\l)dq-^ 

which do commute [E7|, as their classical counterparts. 



Since it is still unknown how to quantize the relation ( |2.29| ) in case of the dynamical 
r-matrices, we again have to rely on good luck trying to find a quantum SoV. An addi- 
tional obstacle is provided by the fact that in our case, even classically, {B{u), B{y)\ 7^ 0. 
Therefore, there is little hope to construct quantum x^ as zeroes of a commuting family 
of operators B{u) like in case of the GL{N) magnet. Instead, we shall rather look for the 
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kernel K{xi,X2^ Qlo'i, Qs) of the integral operator (classically, canonical transformation) 
intertwining the xQ and q representations. 

The first of the classical separated equations (|7.9| ) is easy to quantize. It expresses the 
conservation of the total momentum and allows to eliminate one pair of variables from K: 



Q = 93 



-idr 



{Q - qs)K = 



{Q-q3)K = ^ K = 5{Q - q'i)K{xi,X2\qi,q2,qz), 
[d,, + d,, + d,, + dQ)K = ^ {d,, + d,, + d,,)k = 0, 

k = k{xi, X2|«?13, ?23)- 

To determine the kernel K, let us try to quantize the equations ( |7.7| ). Making the 
substitutions 



id. 



(the last one is a quantum correction found experimentally) one obtains from (|7.7|) the 
system of 4 first order differential equations for K 



(7.11) 



which are easily solved producing the result which most conveniently can be written using 
the variables x± = xi ± X2 and ^± = qis ± 523: 





+ ^913 


)K+{g- 


l)[C(xi) 




- 923) 


-Ciqis) 


+ C(913 


- q23)]K = 







+ ^913 


)k+{g~ 


1)[C(X2) 




- 523) 


-C(913) 


+ C(913 


- q23)]k = 







^923 


)k+{g- 


l)[C(xi) 


-Cixi 


- 913) 


-C(923) 


- C(9l3 


- q23)]k = 







^923 


)k + {g- 


1)[C(X2) 


-ax2 


- 913) 


-C(923) 


- C(9l3 


- q23)]k = 






K = 5{x^- ^+)/C(x+, x_; ^_) 



a 



^-+X. 



cr(xi)cr(x2)o-(^- 



(7.12) 



The above argument, of course, has only heuristic value and provides no guarantee 
that the kernel K thus constructed would produce SoV. What is neccessary to verify is 
that the integral operator with the kernel K transforms an eigenfunction '^{qi, i?2 5 ^3) of the 
hamiltonians if„ ( [7.2| ) satisfying (|2.13| ) into the function \E'(xi, X2, Q) satisfying separated 
equations of the type (|2.14| ). 
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The observation which is crucial for estabhshing SoV is that the kernel K satisfies the 
difi^erential equations 



[-idq - Hl]K = 0, (7.13a) 

-{H; + g{g - iWMx,) - ig{g - l){g - 2)p'(x,))] K = (7.13b) 
where if* is the Lagrange adjoint of Hn 

^{q){H^){q) dq= [ (H* ip){q)i;{q) dq. 



The equations (7.13) can be interpreted as the quantum analog of the equations ( fOj ) 
(note the quantum corrections in g\). Consider now the integral transform 



^{xi, X2,Q) = JJJ dqidq2dqs K{xi, X2, Q\qi, ^2, 93)^(51, ^2, Qs) (7.14) 

Acting on \1/ with the differential operators 

Q = -idQ - /ii, (7.15a) 
Vj = idl^ + hrdl^ - i{h2 + 3^7(^7 - l)p(x,))c>,, 

- (/is + g{g - l)hip{xj) - ig{g - l){g ~ 2)p'{xj)) (7.15b) 

and supposing that '^{q) is an eigenfunction of H^, perform the integration by parts using 
the relations (7.13). The resulting bulk part of the integral is zero. It remains only to 
find such the limits of integration which would not contribute to the result. Omitting the 
details of the guesswork we report only the final result. 

The function \E'(x_|_, x_, Q) obtained from an eigenfunction ^'(^+, q^) of via the 
integral transform 

5(x+,x_,g)= / ^d^_/C(x+,x_;^_)^'(x+,^_,Q) (7.16) 



satisfies the differential equations 

Q$ = 0, Vj^ = 

which imply the SoV 

^{xi,X2,Q) = e'^'^tP{xi)tP{x2] 
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where tl^lx) satisfies a third-order analog of the Lame differential equation 



ii;'" + hitl;" - i{h2 + 3g{g - l)p{x))tl;' 

(7.17) 

- {h3 + gig - l)p{x)hi - ig{g - l){g - 2)p'{x))ij = 0. 

The question of the correct boundary conditions for the separated equation ( [7.17| ) is 
presently under study. 

For the degenerate case of trigonometric potential which corresponds to replacing 
p{q) by 1/ sin^ q in (|7.2| ) and a{u) by sinw in ( [7.4|) the spectrum and eigenfunctions of Hn 
have been well known quite a while ago j^^. The eigenfunctions are labelled by triplets 
z7 = (z^i, z^2, '^s) of integers Uj such that z^i < z^2 < ^^3- The corresponding eigenvalues of 
Hn are given by 

hi = 2{iJi + H2 + Us) (lJ'i = ^i-g 

/i2 = 4(/xi/X2 + + /U2/U3) I 1^2 = ^2 (7.18) 
hs = 8;Ui//2At3 I A*3 = '^s + fi' 

The eigenfunctions \l/i7 have the structure \E',7(^ = \I/ooo(^)^i7(^ where \l/ooo(^) = 
sin^ qi2 sin^ ^13 sin^ ^23 is the vacuum eigenfunction corresponding to u = (0, 0, 0) and 
$/7(g) are symmetric Laurent polynomials in variables tj = e^*'^^ known as Jack polynomials 



The SoV is produced by the same integral kernel K up to replacing a with sin in 
( [7.12| ). The separated eigenfunctions ilj{y) satisfy the differential equation 

+ -^(h2 + 3^^^^") ^' - (hs + ^^^^/ii + 2tg{g-l){g- 2) ^] ^ = 
\ sm X / \ sm X sm x / 

(7.19) 

and can be factorized i/j^ix) = 'ilJoQo{x)fi7{x) into the product of the vacuum factor 
'^ooo{x) = sin^^x and a Laurent polynomial ip^; in variable t = e^*^. Despite the huge 
body of facts known about the Jack polynomials, the last factorization property seems to 
be a new result. 



A more detailed exposition of the above results will be published elsewhere [|9 
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8. Discussion 

The above examples show the diversity of models allowing SoV and give a support to 
the opinion that the domain of SoV method might be very large, even including all the 
models subject to the classical Inverse Scattering Method and their quantum counterparts. 

Let us enlist, in conclusion, some problems whose solution could strengthen the po- 
sitions of SoV. The most obvious object of study is provided by the class of integrable 
models described by numeric unitary /^-matrices, and in the first turn, s/(n)-invariant 
magnets. The s/(2) case being well enough studied, the sl{N) case seems to present only 
calculational difficulties. The case of trigonometric i^-matrices should not differ consid- 
erably from the rational (s/(A^)-invariant) case. For instance, it would be interesting to 
generalise the results of section 5 to the relativistic Toda field theories. In the case of 
elliptic /^-matrices discussed in section 6 the problem of notrivial normalization of the 
Baker- Akhiezer function arises which leads to the complications with the quantization. 
However, the success with the classical XYZ magnet allows to hope that the quantization 
problem would be solved. A similar, but more difficult problem of choosing the correct 
normalization of B-A function arises in case of /^-matrices corresponding to the simple Lie 
algebras other than SL{N). The problem is not yet solved even in the classical case. The 
most difficult problems arise in case of dynamical r-matrices where neither general theory 
exists, no quantization rules are known. 

The problem which may be more important than studying all the particular examples 
is to understand the algebraic structures underlying the SoV and to explain why the "magic 
recipe" of taking the poles of B-A function does work. The recent paper might be the 
first step in that direction. 
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